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1.  Introduction 


The  application  of  force  and  energy  approaches  to  solve  problems  in  mechanics 
have  competed  side  by  side  for  160  years.  Historically,  the  force  balancing  ap¬ 
proach  preceded  the  development  of  energy  methods  that  are  based  on  conservation 
principles.  Because  of  the  relation  that  links  a  change  in  energy  (E)  to  the  action 
of  a  mechanical  force  (F)  acting  through  a  distance,  dE  =  — F  dx ,  we  are  content 
that  the  2  methods  produce  equivalent  results  when  applied  to  the  same  problem. 
The  method  of  choice  depends  on  the  particular  problem  to  be  addressed  and/or 
the  proficiency  of  the  practitioner  with  the  various  methods.  On  one  hand,  the  force 
method  can  be  conceptually  more  straightforward;  however,  an  energy  approach, 
once  properly  conceptualized,  often  proves  more  efficient  in  practice. 

The  tension  between  force  and  energy  approaches,  however,  has  had  an  evolving 
history.  Leaders  in  the  19th  century  physics  community1  debated  long  and  hard 
on  notions  of  conservation,  and  the  relationship  between  force  approaches  and  the 
newer  energy  approaches.  A  prevailing  concept  of  force  “conservation,”  champi¬ 
oned  by  Faraday  and  taken  to  mean  the  “indestructibility  and  transformability  of 
natural  agents,”  was  considered  mathematically  imprecise  and  ambiguous  by  oth¬ 
ers.  Rankine,  on  the  other  hand,  had  posited  the  law  of  energy  conservation  in  1853. 
By  the  end  of  the  decade,  physicists  such  as  Helmholtz,  Thomson,  Rankine,  and 
Maxwell  wished  to  excise  the  force-conservation  concept  from  the  technical  vo¬ 
cabulary. 

In  this  report,  we  examine  1 -dimensional  (1-D)  lattices  from  both  a  force  and  en¬ 
ergy  approach.  The  goal  is  merely  to  confirm  the  compatibility  of  the  2  methods 
as  applied  to  the  mechanical  interactions  of  1-D  lattices.  Specifically,  in  the  energy 
approach,  energy  is  directly  tabulated,  while  force  is  treated  as  a  derivative  quan¬ 
tity,  based  on  the  relation  F  =  —dE/dx.  It  is  determined  that  the  application  of  this 
relation  both  macroscopically  (i.e.,  where  F,  E,  x  are  the  lattice’s  external  force, 
total  energy,  and  body  dimension,  respectively)  and  intemodally  (i.e.,  where  F,  E, 
x  are  the  intemodal  force,  energy,  and  separation,  respectively)  are  compatible  with 
the  force  approach  to  describing  the  lattice  equilibrium. 

The  analysis  of  mass  nodes  interconnected  with  springs  is  a  well-established  field 
and  forms  the  basis  for  the  finite-element  method  (FEM).23  The  use  of  both  force 
and  energy  methods  in  FEM  has  been  employed,  historically.  When  generalizing 
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the  approach  to  large  systems,  the  lattice  connectivity  is  typically  expressed  by 
way  of  a  stiffness  matrix,  to  which  the  methods  of  linear  algebra  are  applied  for  a 
solution.  However,  the  traditional  mechanical  analysis  of  lattices  (and  this  applies  to 
the  FEM  as  well)  is  often  accompanied  by  approximations  of  the  physical  system. 
One  common  simplification  is  the  “harmonic  assumption,”  in  which  the  analysis 
assumes  the  lattice  springs  are  linear  (i.e.,  stiffness  is  independent  of  extension). 
Another  typical  assumption  is  that  of  a  simply  connected  lattice,  in  which  each 
node  of  the  lattice  interacts  with  only  its  nearest  neighbors. 

In  this  report,  however,  we  choose  not  to  limit  ourselves  in  the  manner  described 
above.  This  analysis  considers  fully  connected  1-D  lattices,  in  which  each  node 
of  the  lattice  experiences  a  mechanical  interaction  with  every  other  node  in  the 
lattice.  Thus,  non-local  node  interactions  are  considered.  In  addition,  no  harmonic 
assumption  is  made,  in  which  the  spring  constants  are  invariant  with  extension.  In 
fact,  no  particular  spring-force  relationship  is  assumed  at  all,  except  for  the  purposes 
of  illustration.  Our  only  assumption,  for  the  purpose  of  these  derivations,  is  that 
the  component  spring  force  f(x)  is  knowable,  from  which  the  necessary  external 
force(s)  may  be  determined  to  produce  a  given  value  of  nominal  lattice  spacing  A. 

In  the  iV-node  1-D  lattices  that  are  considered  in  this  report,  discussion  is  limited 
to  symmetric  equilibrium  configurations.  If  the  notation  Fj  is  used  to  denote  the 
Cartesian  external  load  applied  to  node  j  of  the  lattice  (positive  to  the  right),  then 
the  symmetric  equilibrium  condition  is  simply  that  Fj  =  —  F^+i-j- 

2.  The  Force  Approach 

The  force  approach  is  typically  one  in  which  known  external  forces  are  applied  to  a 
body,  which,  in  turn,  establishes  resultant  internal  forces  in  the  body.  If  the  consti¬ 
tutive  laws  governing  the  body  are  known  (which  in  the  case  of  a  lattice  means  the 
spring  force  relation  f(x)),  then  a  deformation  in  the  body  is  likewise  established, 
associated  with  the  applied  external  load.  In  the  current  analysis,  however,  I  use  the 
known  constitutive  behavior  of  the  lattice  to  instead  answer  the  question  of  what 
force  must  be  externally  applied  in  order  to  establish  a  given  deformation. 

The  method  for  calculating  the  internal  forces  is  the  ‘Tree-body  diagram,”  as  de¬ 
scribed  in  any  undergraduate  book  on  physics4  or  statics.5  In  the  case  of  a  non¬ 
accelerating  system,  the  body  is  conceptually  “cut  in  two”  and  the  internal  forces 
across  the  cut  of  the  body  must  be  of  a  magnitude  and  direction  to  precisely  bal- 
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ance  the  external  forces  being  applied  to  that  section  of  the  body.  In  general,  this 
force  balancing  requires  vectorial  addition;  however,  because  the  problem  under 
consideration  is  a  1-D  lattice,  the  force  balance  becomes  a  purely  scalar  process. 


2.1  The  2-Node  Structure 

The  archetypal  2-node  interaction  is  shown  in  Fig.  1.  In  this  interaction,  the  node 
pair  1tA2  has  been  brought  together  to  a  separation  distance  of  A,  requiring  the  ap¬ 
plication  of  symmetric  forces  F\  and  F>.  The  function  f(x).  here  and  throughout 
this  report,  is  taken  as  the  internal  (compressive)  spring  force  that  acts  between 
any  2  arbitrary  nodes  of  the  lattice,  when  separated  by  a  distance  x  (note:  x  is  not 
the  displacement  from  a  rest  configuration,  but  rather  the  physical  spring  length!). 
While  the  convention  is  here  employed  that  a  positive  force  is  compressive,  tensile 
situations  are  covered  with  no  loss  in  generality. 


— — , 


A 


Fig.  1  A  1-D,  2-node  structure  with  spacing  A 

In  this  figure  and  throughout  this  report,  I  use  the  notation  that  fk  =  f(k A)  rep¬ 
resents  the  “spring”  force  between  a  pair  of  nodes  that  are  separated  by  a  distance 
of  sk  =  kA  (in  general,  k  need  not  be  an  integer).  Thus,  in  Fig.  1 ,  f\  =  /(A). 
Symmetric  equilibrium  dictates  that  F\  =  —F2. 

In  this  simplest  of  interactions,  conceptually  cutting  the  lattice  between  nodes  1  and 
2  would  reveal  that  the  internal  force  across  the  cut,  call  it  F,  is  comprised  solely  of 
the  single  spring  force  f\.  Enforcing  a  force  balance  over  the  cut  “free  body”  (such 
that  the  external  force  F\  must  exactly  balance  the  internal  force  F)  simply  reveals 
that 

Fi  =  F  =  /i  .  (1) 

Thus,  the  separation  distance  x  will  adjust  itself  under  the  symmetrically  applied 
load,  so  as  to  precisely  generate  a  spring  force  of  /i  =  F\  when  x  =  A.  The  internal 
force  F  is  uniform  throughout  the  extent  of  the  2-node  structure.  Because  we  make 
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no  assumptions  about  the  functional  nature  of  the  spring  force  f(x),  no  more  can 
be  said  about  the  relation  between  F  and  A. 


2.2  The  3-Node  Lattice 

The  3-node  structure  depicted  in  Fig.  2  is  the  simplest  geometry  that  can  be  called 
a  “lattice,”  because  of  the  repetition  of  cell  geometry.  Symmetrical  equilibrium  is 
enforced  as  a  constraint,  such  that  F2  (not  shown  in  figure)  is  exactly  zero,  and  F\  = 
— F3.  Because  of  topologically  symmetric  stiffness  about  node  2,  the  uniformity  of 
spacing  between  node  pairs  1  v--2  and  2 -a  3  is  assured,  denoted  here  as  A. 


Fig.  2  A  fully  connected,  1-D,  3-node  lattice  with  uniform  spacing  A 

Accounted  for  here  is  the  effect  of  non-nearest-neighbor  interactions,  in  which  node 
pair  1  aA3  is  likewise  joined  by  way  of  a  “spring”  interaction — the  same  f(x)  spring 
that  governs  the  interactions  of  node  pairs  1tA2  and  2vA.  However,  in  the  case  of 
node  pair  1tA3,  the  spring  force  /  is  evaluated  at  a  separation  of  2A,  thus  /2. 

Whether  one  conceptually  cuts  the  structure  between  node  pair  1tA2  or  2vl3.  the 
internal  force  is  found  to  be  /j  +  /2.  Performing  the  balance  between  internal  and 
external  force  leads  to  the  result 

Fi  =  F  =  f1  +  f2  .  (2) 

A  quick  comparison  to  Eq.  1  reveals  that  the  external  force  Fx,  corresponding  to 
a  given  spacing  A  for  the  2-  and  3-node  cases,  is  different,  even  when  employing 
the  same  f(x)  spring  in  both  cases.  Such  a  discrepancy  will  manifest  when  spring 
connectivity  includes  non-nearest  neighbors. 
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2.3  Uniformly  Spaced  Lattice  vs.  Simply  Loaded  Lattice 

To  extend  this  process  beyond  the  3-node  lattice,  a  choice  needs  to  be  made,  for 
a  basic  incompatibility  arises.  On  one  hand,  if  the  lattice  remains  simply  loaded 
at  the  lattice  boundary,  the  nodal  spacing  can  no  longer  remain  uniform.  On  the 
other  hand,  enforcing  a  uniform  lattice  spacing  implies  that  some  form  of  external 
loading  must  be  applied  even  to  internal  nodes  (one  need  not  concern  oneself  with 
how  that  might  be  brought  to  bear). 

Both  possibilities  are  treated  in  the  following  subsections.  Considered  are  uniformly 
spaced  4-,  5-,  and  6-node  lattices,  as  well  as  simply  loaded  4-  and  5-node  lattices. 
However,  some  new  nomenclature  is  required  to  describe  the  added  complexity. 

In  the  case  of  a  uniformly  spaced  lattice,  external  forces  are  necessarily  applied  to 
internal  nodes.  With  the  external  load  being  distributed  over  several  nodes,  the  net 
force  being  applied  to  the  lattice  may  be  defined  as 

N/2 

Fnet  —  y ^Fj  .  (3) 

3= 1 

Note  that,  when  N  is  odd,  it  does  not  matter  whether  the  sum  index  j  is  rounded  up 
or  down  since,  because  of  symmetric  equilibrium,  the  central  lattice  node  will  have 
0  external  force  being  applied  to  it. 

Because  external  forces  are  applied  to  internal  nodes  for  the  case  of  a  uniformly 
spaced  lattice,  no  longer  is  F  sufficient  to  completely  characterize  the  internal  lat¬ 
tice  force.  Thus,  Fj  is  used  to  denote  the  internal  force  in  a  lattice  that  exists  between 
the  node  pair  j  v-  (j  +  1)  (e.g.,  F\  is  the  total  internal  force  that  exists  between  node 
pair  1tA2).  Note  that,  because  of  symmetric  equilibrium,  Fj  =  FN-r 

The  nomenclature  Sk  has  already  been  introduced  as  shorthand  for  the  distance 
k\.  In  a  related  way,  the  syntax  of  rr-x  is  provided  to  mean  the  physical  separation 
between  node  pair  i^rri  of  the  lattice.  For  the  special  case  of  a  uniformly  spaced 
lattice,  rr-x  =  S|m_j|  =  \m—i\X.  Whereas  the  notation  provides  the  specificity  of 
both  nodes  in  the  node  pair,  the  notation  Sk  allows  the  convenience  of  tracking  only 
the  separation  distance  in  node  pair  i/m,  without  the  need  to  know  the  specific 
values  of  i  and  m. 

The  presence  of  a  spatially  dependent  internal  force  allows  for  the  consideration  of 
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an  average  internal  force  within  the  lattice.  In  the  most  general  consideration, 


However,  when  the  lattice  is  uniformly  spaced,  the  above  form  is  simplified,  by 
knowing  that  j+jx  =  X: 


2  _ _ * 

— -  ^  Fj  ,  (uniform  lattice  spacing). 


In  the  case  of  a  simply  loaded  lattice,  however,  the  spacing  within  the  lattice  will 
not  remain  uniform.  Thus,  we  introduce  the  nomenclature  of  A  as  a  fraction  of  A 
representing  the  deviation  in  lattice  spacing  from  the  nominal  value.  In  this  way, 
/i_ a  represents  the  spring  force  evaluated  at  a  nodal  separation  of  A(1  —  A).  With  a 
suitably  large  lattice,  one  would  require  indexed  A j  to  represent  different  deviations 
from  the  nominal  lattice  spacing  for  different  node  pairs.  However,  for  the  scope  of 
this  report,  a  single,  unindexed  A  is  sufficient. 

2.4  The  Uniformly  Spaced  4-Node  Lattice 

In  Fig.  3  is  found  a  4-node  lattice  in  which  uniform  spacing  is  enforced.  In  order  to 
maintain  uniform  spacing,  there  must  be,  in  addition  to  the  forces  F1  and  F4  applied 
at  the  external  boundaries,  forces  F2  and  F:>  applied  upon  the  internal  nodes. 


Conceptually  cutting  the  lattice  between  each  successive 

forces  in  the  lattice  are  revealed  as 

node  pair,  the  internal 

A  =  /i+  h  +  h  , 

(6) 

F2  =  fi  +  2f2  +  f3  , 

(7) 

and 

Kr3> 

CO  J 

(8) 

6 


Fig.  3  A  fully  connected,  1-D,  4-node  lattice  with  uniform  spacing  A 


With  this  knowledge,  the  average  internal  lattice  force  may  be  determined  as 

F  —  -  (^Fi  +  F2  +  F3 
—  fl  +  g/2  +  fs 

If  one  performs  the  internal/extemal  force  balance  with  each  instance  of  lattice 
“cut,”  the  free-body  diagrams  reveal  that 


F\  —  Fi  ,  F2  —  Fi  +  F2  ■  (10) 

Solving  for  the  external  forces  reveals 

Fi  =  Fi  =/i  +  /2  +  /s  ,  (11) 

F2  =  F2-Fi=  f2  (12) 
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with  F4  =  —F\  and  F:i  =  —  F2.  The  net  externally  applied  load  is 


Fnet  ~  Fi  +  F2  —  F2 

=  fl  +  2/2  +  h 


(13) 


Because  the  net  external  load  is  not  uniformly  present  across  the  lattice,  it  is  seen 
to  differ  from  the  average  internal  load,  F,  given  in  Eq.  9. 


As  a  purely  illustrative  example,  consider  a  spring  with  the  relation  f(x)  oc  1/a;2. 
For  such  a  hypothetical  example,  the  significance  of  F2  compared  to  the  net  load 
Fnet  can  be  evaluated  from  Eqs.  12  and  13  as 


F2  _  1/22 

Fnet  ~  l/l2  +  2/22  +  1/32 


9 

58 


«  0.156 


(14) 


For  this  hypothetical  example,  the  significance  of  F2  relative  to  Fnet  is  indepen¬ 
dent  of  A.  That  is  by  no  means  guaranteed,  when  considering  cases  involving  more 
complex  functional  relationships  for  f(x). 


2.5  The  Uniformly  Spaced  5-Node  Lattice 

Consider  the  uniformly  spaced  5-node  lattice  in  Fig.  4.  It  may  be  useful  to  note  that 
the  total  number  of  springs  in  a  fully  connected  Ar-nodc  lattice  is  Ar(Ar  —  l)/2; 
namely,  N  —  1  springs  with  a  separation  of  s4,  N  —  2  springs  with  a  separation  of 
s2,  ...  ,  and  1  spring  with  a  separation  of  S(jv-i). 


As  before,  conceptual  cuts  are  made  through  the  lattice  between  each  successive 
node  pair,  in  order  to  reveal  the  internal  loads: 


Fi  =  f  1+  /2+  fs  +  f 4  , 

F2  —  f\  +  2/2  +  2/3  +  /4  , 

and 

F3  =  F2  ,  F4  =  Fi  . 

From  these,  the  average  internal  force  may  be  calculated  as 

F  =  |  (A  +  F2  +  F3  +  F4) 

—  fi  +  2^2  +  2^3 


(15) 

(16) 


(17) 


(18) 
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A  A  A  A 


Fig.  4  A  fully  connected,  1-D,  5-node  lattice  with  uniform  spacing  A 

A  force  balance  at  each  instance  of  lattice  cut  reveals  the  internal  forces: 

Fi  —  Fi  ,  F2  =  Fi  +  F2  . 

These  2  equations  may  be  solved  for  the  external  forces: 

F1  =  F1  =fi  +  f2  +  h  +  h  , 

F2  =  F2  —  Fi  =  f2  +  f3  , 

with  F5  =  —F\ ,  F4  =  —  F2  and  F:>  =  0,  by  way  of  symmetric  equilibrium 


(19) 


(20) 

(21) 

.  The  net 
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externally  applied  load  is,  therefore, 


Fnet  —  Fi  +  F2  —  F2 

—  f\  +  2/2  +  2/3  +  /4 


(22) 


The  net  applied  lattice  force,  Fnet,  is  seen  to  differ  from  the  average  internal  force, 
F,  in  the  lattice,  given  by  Eq.  18. 


As  before,  we  can  compute  the  significance  of  F2  relative  to  Fnet,  employing  a 
hypothetical  spring  for  which  f(x)  oc  1/x2.  The  result  is 


F-2  _  1/22  +  1/32 

Ket  ““  l/l2  +  2/ 22  +  2/ 32  +  1/42 


52 

257 


«  0.202 


(23) 


2.6  The  Uniformly  Spaced  6-Node  Lattice 

The  final  case  of  a  uniformly  spaced  lattice  considered  in  this  report  is  that  com¬ 
prised  of  6  nodes,  as  depicted  in  Fig.  5.  Making  conceptual  cuts  through  the  lattice 
between  each  successive  node  pair  is  done  to  evaluate  the  internal  loads: 


Fi  =  h+  h+  h+  U  +  f 5 

(24) 

-^2  =  fi  +  2/2  +  2/3  +  2/4  +  /5 

(25) 

-^3  =  fi  +  2/2  +  3/3  +  2/4  +  f5 

(26) 

h  =  F2  ,  F5  =  F1  . 

(27) 

As  in  prior  configurations,  the  average  internal  force  may  be  evaluated  as 


F  —  -  (^F\  +  F2  +  F3  +  F4  +  F$j 
8  9  8 

—  fi  +  7/2  +  7/3  +  7/4  +  fh 


(28) 


A  force  balance,  performed  at  each  conceptual  lattice  cut,  reveals  the  magnitude  of 
the  internal  forces: 


F\  —  F\  ,  F2  —  F\  +  F2,  ,  T3  —  F\  +  F2  +  F3  .  (29) 
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f3 


f3 

— 's/v 

u  I 


2 


Fig.  5  A  fully  connected,  1-D,  6-node  lattice  with  uniform  spacing  A 
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These  3  equations  may  be  solved  for  the  applied  external  forces: 


Fl=Fl  =/l  +  /2  +  /3  +  /4  +  /5  ,  (30) 

F2  =  F2-F1  =  /2  +  /3  +  /4  ,  (31) 

F3  =  F3-F2  =  f3  ,  (32) 


with  Fq  =  —Fi,  F5  =  —F2,  and  F4  =  —F3.  The  net  external  load  (again  distinct 
from  F)  is  therefore 


Fnet  —  F\  +  F2  +  F3  —  F3 

—  fl  +  2/2  +  3/3  +  2/4  +  /5 


(33) 


For  a  hypothetical  lattice  spring  in  which  f(x)  oc  1/a;2,  the  relative  significance  of 
both  F2  and  F3  may  be  evaluated: 


F2  _  1/22  +  1/32  +  1/42  _  1525 

~  l/l2  +  2/22  +  3/32  +  2/42  +  1/52  “  7194 


(34) 


and 


= _ 1/3^ _ 

Fnet  l/l2  +  2/22  +  3/32  +  2/42  +  1/52 


400 

7194 


«  0.056 


(35) 


Of  course,  a  different  functional  relationship  for  the  lattice  spring  would  imply  dif¬ 
ferent  quantitative  values  for  the  significance  of  F2  and  F3. 


2.7  The  Simply  Loaded  4-Node  Lattice 

Consider  now  a  lattice  that  is  simply  loaded  at  its  boundaries.  While  the  internal 
load  Fj  may  be  evaluated  by  conceptually  cutting  the  lattice  between  any  2  adjacent 
nodes  j  and  j  +  1,  the  absence  of  external  forces  at  the  internal  nodes  implies  that 
the  internal  load  is  everywhere  uniform,  or  Fj  =  F.  Figure  6  depicts  such  a  lattice, 
comprising  4  nodes.  Because  the  internodal  force  relationships  differ  for  the  internal 
nodes  (2,  3)  vis-a-vis  the  boundary  nodes  (1,4),  the  interlattice  spacing  of  node  pairs 
1  tA2  and  3^4  will  differ  by  an  amount  A  A  from  that  of  node  pair  2vA.  The  lattice 
spacing  of  node  pair  2^3  is  designated  A,  since  that  pair  is  the  farthest  from  the 
influence  of  boundary  effects. 
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(1-AU  A.  (1  —A)  A 


Fig.  6  A  fully  connected,  1-D,  4-node  lattice  with  nominal  spacing  A  and  simple  boundary 
loading 


Performing  a  conceptual  cut  between  each  successive  node  pair  so  as  to  tabulate  the 
internal  forces  gives 

F\  —  fi—A  +  /2-A  +  /3-2A  ,  (36) 

F-2  —  /1  +  2/2_a  +  /3-2A  ,  (37) 


with  F3  =  F],  However,  in  the  case  of  a  simply  loaded  lattice,  the  internal  loads 
must  be  everywhere  uniform.  Thus,  F  =  F\  =  F2,  leading  to  the  constraint  on  A 
that 

,/l  +  h-  A  =  fl—A  ■  (38) 

The  external  load  F1?  being  equal  to  F,  can  therefore  be  expressed  in  several  ways: 


Fi  =  F 


/l— A  +  /2-A  +  /3  — 2A 

fl  +  2/2_a  +  /3-2A 


(39) 


Only  if  the  spring  function  is  explicitly  defined  may  the  value  of  A  be  quantita¬ 
tively  ascertained.  If,  for  example,  the  lattice  spring  obeyed  f(x)  oc  1/x2.  then  the 
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constraint,  Eq.  38,  would  take  the  form 


1  1  _  1 

l2  +  (2- A)2  “  (1- A)2  ’ 


(40) 


which  may  be  numerically  solved  for  A  «  0.117  (limiting  A  to  real  values  in  the 
range  A  <  1).  Thus,  an  11.7%  variation  in  lattice  spacing  would  be  evident  across 
a  simply  loaded  4-node  lattice,  for  which  f(x)  oc  l/x1  governed  the  spring  force. 


2.8  The  Simply  Loaded  5-Node  Lattice 


The  final  lattice  that  is  considered  in  this  report  is  a  simply  loaded  5-node  lattice 
(Fig.  7). 


(1  —A)  A  /l  A  (1-AU 


Fig.  7  A  fully  connected,  1-D,  5-node  lattice  with  nominal  spacing  A  and  simple  boundary 
loading 


The  force  balance  may  proceed  at  each  conceptual  cut  across  the  lattice,  between 


14 


successive  nodes: 


A  —  fl-A  +  /2-A  +  /3-A  +  /4-2A  ,  (41) 

A  =  /1  +  /2  +  /2-A  +  2/3_a  +  /4-2A  ,  (42) 


and 

A  =  A  j  A  =  A  •  (43) 

Since  a  uniform  internal  load  is  manifested  throughout  the  lattice,  it  is  the  case  that 
F  =  Fi  =  F2,  which  leads  to  the  constraint  on  A  that 


fl—A  ~  fl  +  fl  +  /3-Z 


(44) 


The  external  load  F1,  being  equal  to  F,  can  be  expressed  in  several  ways: 


Fi  —  F  — 


fl-A  +  /2-A  +  /3-A  +  /4-2A 

fl  +  /2  +  /2-A  +  2/3_a  +  /4-2A 


(45) 


If  a  spring  function  is  specified,  the  value  of  A  will  be  determinable.  Considering  a 
hypothetical  spring  for  which  /(x)  oc  1/x2,  the  constraint,  Eq.  44,  becomes 


(1  —  A)2  l2  22  (3  -A)2  ’ 


(46) 


which  can  be  numerically  evaluated  to  ascertain  that  A  «  0.147  (limited  to  real 
solutions  for  which  A  <  1). 


3.  The  Energy  Approach 

In  Section  2,  a  force  approach  was  used  to  develop  the  equations  of  equilibrium  in  a 
series  of  1-D  lattices.  When  nodes  in  the  lattice  are  not  accelerating,  equilibrium  is 
established,  with  the  help  of  free-body  diagrams,  by  developing  a  balance  between 
internal  and  external  forces.  A  general  formulaic  tabulation  of  equilibrium  for  a  1-D 
lattice  of  arbitrary  size  was  not  sought  using  the  force  approach.  Rather,  individual 
1-D  lattices,  ranging  from  2  to  6  nodes,  were  examined,  under  2  separate  behavioral 
assumptions:  that  of  constant  lattice  spacing  and  that  of  simple  loading  at  the  lattice 
boundary.  The  process  was  a  tedious  one,  in  which  a  full  accounting  of  internodal 
interactions  was  ascertained  and  tabulated.  Were  the  lattices  in  2  or  3  dimensions 
rather  than  1 ,  the  formulations  would  be  still  more  intricate,  as  vectorial  calculations 
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would  be  required  for  component  resolution. 


In  the  force  approach,  the  concept  of  energy  was  totally  absent  from  the  analysis.  In 
this  section,  the  goal  is  to  re-examine  those  same  lattices  using  an  energy  approach, 
in  which  force  is  a  quantity  derived  from  the  energy.  The  goal  is  to  establish  a 
compatible  result. 


In  a  lattice,  let  the  potential  energy  that  any  2  nodes  possess  as  a  result  of  their 
mutual  interaction  at  a  separation  distance  of  x  be  given  as  e(x).  If  the  nodes  are 
materially  indistinguishable,  then  half  of  that  potential  energy  may  be  partitioned 
to  each  of  the  2  interacting  nodes.  The  energy  associated  with  node  i  of  an  iV-node 
lattice  is,  therefore,  the  sum  of  its  interactions  with  all  the  other  nodes  of  the  lattice: 

1  N 

E,  =  j  £  e(»  ,  (47) 

m=  1 


where  e(0)  is  taken  as  0  to  address  the  interaction  of  node  i  with  itself,  and  where 
the  factor  of  1/2  accounts  for  the  partition  of  potential  energy  to  both  nodes  partici¬ 
pating  in  an  interaction  ( i.e .,  only  1/2  the  interaction  energy  belongs  to  node  i).  The 
total  energy  of  the  lattice  is  obtained  by  repeating  this  process  over  all  the  nodes  of 
the  lattice. 

1  N  N 

e=~2T,T, e<»  •  <48> 

i= 1  m= 1 

The  macroscopic  force  associated  with  changing  the  dimension  of  the  overall  lattice 
structure  is  derived  from  the  total  energy  as 


dE 

dL 


(49) 


whereas  the  force  of  single  internodal  interaction  at  a  spacing  of  s  is,  correspond¬ 
ingly, 


de 

ds 


(50) 


The  key  question  examined  in  this  report  is  determining  whether  these  derived 
forces  (Eqs.  49  and  50)  in  the  context  of  Eq.  48  are  compatible  with  the  earlier 
derived  force  approach  of  Section  2. 


A  prerequisite  in  addressing  this  question  is  understanding  the  relationship  between 
a  change  in  the  global  lattice  size  (L)  and  the  corresponding  change  in  the  lattice 
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spacing  of  node  pair  iv'-m,  which  is  given  by  ™x.  To  do  so,  the  chain  rule  is  applied. 


dei^x)  de(™x)  d™x  dX 

dL  d™x  dX  dL 

The  first  term  of  the  right-hand  side  is  simply 

dei^x) 


(51) 


dEx 


=  -f  C»  , 


(52) 


by  definition.  To  the  first  order,  the  last  term  of  Eq.  51,  under  simple  expansion,  is 


dX  _  X 
dL  ~  L 


(53) 


which  is  independent  of  the  node  pair  is-rn.  Equation  51  may  thus  be  employed 
when  taking  the  derivative  of  the  energy  equation,  Eq.  48,  with  respect  to  the  body 
dimension  L,  to  obtain  the  macroscopic  force: 


F 


X 

2 L 


N  N 


£  £  /(’» 


dnlx 

~dX 


(54) 


Similar  to  energy,  /( 0)  =  0  to  remove  any  effect  of  self  force.  The  multiplication  of 
the  local  force  by  dr,-  x/dX  is  key  to  a  proper  force  enumeration  and  has  been  seen 
before  in  my  previous  work  on  lattice  force6-8  (as  ds/dX). 


3.1  Nodal  Forces 

Whereas  Section  3  provides  the  means  to  evaluate  the  global  lattice  force  in  terms 
of  the  spring-component  forces,  by  way  of  an  energy  derivative,  it  does  not  provide 
the  means  to  evaluate  the  component  nodal  forces.  To  obtain  the  externally  applied 
nodal  forces  from  the  energy  approach,  one  has  to  imagine  a  connected  lattice  in 
which  1  of  the  nodes,  i,  is  displaced  from  its  equilibrium  position  by  an  infinites¬ 
imal  amount,  Ade*.  The  pre-existing  applied  nodal  force  that  moves  through  that 
displacement  causes  an  energy  change  such  that 


Fi 


dE 
X  det 


(55) 


Since  E  is  composed  of  a  sum  of  internodal  e  terms,  one  can  evaluate  the  derivative 
using  a  chain  rule: 


de(™x)  de(™x)  d™x 

X  dei  d™x  X  det 


(56) 


17 


Therefore, 


N  N 


X 


d™x 
A  dei 


(57) 


t=l  m=  1 

When  neither  of  the  indices,  i  nor  m,  correspond  to  the  displaced  node  i,  the  dis¬ 
tance  ™x  remains  unchanged  and  thus,  d™x/dei  —  0  (i  ^  i,  m).  Also,  the  distance 
n-x  does  not  depend  on  the  order  of  i  and  m,  so  that  ™x  =  ^x.  This  symmetry,  and 
the  requirement  that  either  i  or  m  take  on  the  value  of  i  for  a  non-zero  term  to  arise, 
allows  for  a  condensation  of  Eq.  57: 


N 

E 

m= 1 


f,  =  -  > :  /<» 


d™x 

A  dej 


(58) 


Note  that,  because  e  is  infinitesimal,  one  may  safely  approximate  fk±e  as  ft  in 
the  evaluation  of  /(r"x)  terms  in  Eq.  58.  However,  given  an  infinitesimal  positive 
displacement  of  node  i,  the  final  term  of  Eq.  58  becomes 


d™;x 
A  det 


{1  (m  <  i ) 
0  (m  =  i ) 
—  1  (m  >  i ) 


(59) 


When  a  lattice  is  uniformly  spaced,  the  resultant  Ft  indicates  the  external  force 
being  applied  to  node  i  in  order  to  maintain  the  uniform  spacing. 

For  the  case  of  the  simply  loaded  lattice,  the  boundary  load  may  be  determined  by 
choosing  i  —  1  (or  i  —  N)  for  the  evaluation  of  Eq.  58.  However,  when  i  is  chosen 
as  an  internal  node  of  a  simply  loaded  lattice,  it  is  known  that  the  applied  force  on 
internal  nodes  is  zero.  Therefore,  setting  the  resultant  force  equivalent  to  0  instead 
generates  a  constraint  equation  on  the  internodal  forces,  of  the  kind  indicated  in 
Eqs.  38  and  44.  These  constraints  can  be  used,  once  a  lattice- spring  force  relation  is 
specified,  to  evaluate  the  nonuniformity  (A)  of  the  equilibrium  lattice  spacing  for  a 
simply  loaded  lattice. 

We  now  proceed  to  reanalyze  the  lattices  of  Section  2,  in  light  of  Eqs.  54  and  58. 
For  the  analysis,  the  shorthand  ek  =  e(sk)  =  e(kX)  is  introduced,  in  a  fashion 
analogous  to  fk. 
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3.2  Simplification  for  Uniformly  Spaced  Lattice  Structures 


For  a  uniformly  spaced  lattice,  in  which  n-x  =  s\m-i\  =  \m  —  i\\,  the  second  and 
last  terms  in  the  derivative  chain  of  Eq.  51  immediately  follows  as 


and 


d  i x  ds\m_i\ 

d\  d\ 


m  —  i  |  ,  (uniform  lattice  spacing) 


(60) 


dX  X  1 
dL  =  L  =  N-  1 


,  (uniform  lattice  spacing). 


(61) 


On  the  other  hand,  when  the  lattice  is  simply  loaded,  Eq.  60  needs  to  be  modified 
to  account  for  the  A  terms  associated  with  the  nonuniformity  of  lattice  spacing  for 
the  node  pair  Am.  Nonetheless,  under  simple  expansion,  to  the  first  order 


dmiX  _  ™x 

~dX  ~  T 


(62) 


Likewise,  Eq.  61  will  necessarily  revert  to  Eq.  53  for  a  simply  loaded  lattice — 
nonetheless,  X/L  remains  a  parameter  independent  of  i  and  m,  accounting  for  its 
removal  from  the  double  summation  of  Eq.  54. 


3.3  The  2-Node  Structure 

For  this  discussion,  please  refer  to  Fig.  1 .  The  simplifications  of  Section  3.2  apply  to 
this  case.  As  the  lattice  size  ( N )  increases  (in  later  report  sections),  it  helps  organize 
the  mind  first  by  tabulating  the  nodal  energies  (Eq.  47): 


2  E\  —  e\ 

2E2  =  e\ 

2 

E  =  YJEl  =  el  (63) 


A  direct  application  of  Eq.  54  gives 


F  = 


=  fi 


(64) 


which  matches  Eq.  1  as  expected. 
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3.4  The  3-Node  Lattice 


For  this  discussion,  please  refer  to  Fig.  2.  The  simplifications  of  Section  3.2  apply 
to  this  case.  The  nodal  energy  tabulation  (Eq.  47)  is 


2E\  —  e\  +  e2 
2E2  =  2ei 
2E3  =  e\  +  e2 

3 

E  =  'y  ^  Ei  =  2ei  +  e2  (65) 


Recall,  the  subscript  on  E%  refers  to  a  node  number,  whereas  the  subscript  on  ep- 
refers  to  internodal  separation  as  a  multiple  of  A.  A  direct  application  of  Eq.  54 
gives 


F 


2/i  •  1  +  f2  ■  2 


(66) 


=  fl  +  f2 


which  matches  Eq.  2,  as  derived  using  a  force  approach. 


3.5  The  Uniformly  Spaced  4-Node  Lattice 

For  this  discussion,  please  refer  to  Fig.  3.  The  simplifications  of  Section  3.2  apply 
to  this  case.  A  tabulation  of  the  nodal  energies  (Eq.  47)  gives 


2E\  =  ei  +  e2  +  e3 

2E2  =  2ei  +  e2 

2E%  =  2ei  +  e2 

2T/4  =  ei  +  e2  +  e3 
4 

E  =  y  ^  Ej  =  3ei  +  2e2  + 
%= 1 

A  direct  application  of  Eq.  54  gives 


F 


2  |^3/i  •  1  +  2/2  •  2  +  /3  •  3 

fi  +  3/2  +  /3 


(67) 


(68) 


A  comparison  to  the  force  approach  of  Section  2.4  gives  an  interesting  revelation. 
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Because  a  uniformly  spaced  lattice  has  multiple  points  of  external  force  application, 
there  is  not  a  simple  value  of  F  against  which  to  compare  Eq.  68.  Instead,  there  is 
the  net  applied  force  (Fnet)  of  Eq.  13  and  the  average  internal  force  (F)  of  Eq.  9. 
Interestingly,  this  value  of  macroscopic  force  F  given  by  Eq.  68  matches  F.  the 
average  internal  force  inside  the  lattice. 


To  generate  the  applied  nodal  forces,  the  relation  Eq.  58  needs  to  be  applied  to 
successive  nodes  of  the  uniformly  spaced  4-node  lattice.  A  direct  application  to 
node  1  gives 


F\  = 


/i  •  (-i)  + /2  •  (-1)  + /3  •  (-i) 


(69) 


which  simplifies  to 


F\  =  ./,  +  /2  +  ,/3  , 


(70) 


which  matches  the  expression  given  in  Eq.  1 1  for  the  node  1  loading  of  the  uni¬ 
formly  spaced  4-node  lattice.  A  direct  application  to  node  2  gives 


F2 


h  ■  (1)  +  fi  •  (-1)  +  f-2  •  (-1) 


=  /2 


(71) 


which  matches  the  expression  given  in  Eq.  12  for  the  node  2  loading  of  the  uni¬ 
formly  spaced  4-node  lattice. 

3.6  The  Uniformly  Spaced  5-Node  Lattice 

For  this  discussion,  please  refer  to  Fig.  4.  The  simplifications  of  Section  3.2  apply 
to  this  case.  A  tabulation  of  the  nodal  energies  (Eq.  47)  gives 

2  E\  =  e\  +  e2  +  +  e\ 

2E2  =  2e  i  +  e2  + 

2E%  =  2e  i  +  2e2 
2E/4  =  2e  1  +  e2  +  e% 

2E§  =  ei  +  62  +  63  +  64 

5 

E  =  'y  ^  E{  =  4ei  +  3e2  +  2e3  +  64  (72) 

i= 1 

The  pattern  in  the  term  multipliers  can  be  inferred  as  the  number  of  springs 
present  at  a  distance  of  1A,  2 A,  3A,  etc.,  respectively,  in  the  given  lattice  configura- 
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tion.  A  direct  application  of  Eq.  54  gives 


F 


-  [4/i  •  1  +  3/2  •  2  +  2/3  •  3  +  fA  ■  4 

fl  +  2^2  +  2^3  +  fi 


(73) 


As  in  the  case  of  the  4-node  lattice,  the  value  of  macroscopic  force  F  given  by 
Eq.  73  matches  F,  the  average  internal  force  inside  the  5-node  lattice,  as  given  by 
Eq.  18. 


To  generate  the  applied  nodal  forces,  the  relation  Eq.  58  needs  to  be  applied  to 
successive  nodes  of  the  uniformly  spaced  5-node  lattice.  A  direct  application  to 
nodes  1-3,  respectively,  gives 


Fi 


/i-(-l)  +  /2-(-l)  +  /3-(-l)  +  /4-(-l) 


F2 

F3 


fl  +  /2  +  /3  +  /4 

-  /1  ■(!)  +  /!  •(-!)  +  /2  •(-!)  +  /3  -  (-1) 

/2  +  /3 

-  /2  •  (1)  +  fi  ■  (1)  +  fi  •  (-1)  +  /2  •  (-1)' 


=  0 


(74) 

(75) 

(76) 


matching  the  expressions  given  in  Eqs.  20-21  for  the  loading  of  the  uniformly 
spaced  5-node  lattice.  An  application  to  node  3  gives  the  expected  result  that  F3  = 
0. 

3.7  The  Uniformly  Spaced  6-Node  Lattice 

For  this  discussion,  please  refer  to  Fig.  5.  The  simplifications  of  Section  3.2  apply 
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to  this  case.  A  tabulation  of  the  nodal  energies  (Eq.  47)  gives 


2Fi  = 

=  e-i 

+ 

e2 

+ 

e3 

+ 

e4  +  es 

2E2  = 

=  2ei 

+ 

e2 

+ 

e3 

+ 

e4 

2F3  = 

=  2ei 

+  2e2 

+ 

e3 

2  F4  = 

=  2ei 

+  2e2 

+ 

e3 

2  F5  = 

=  2ei 

+ 

e2 

+ 

e3 

+ 

e4 

2Fg  = 

=  ei 

+ 

e2 

+ 

e3 

+ 

e4  +  es 

E  = 

6 

Ei 

= 

5e 

1  4 

■  4e2  +  3e3  +  2e4  +  e§ 

i= 1 


A  direct  application  of  Eq.  54  gives 


F 


g  |5/i  •  1  +  4/2  ■  2  +  3/3  •  3  +  2/4  •  4  +  /5  •  5 
8  9  8 

/l  +  7_/2  +  7/3  +  r  /4  +  /5 


(77) 


(78) 


As  in  the  case  of  the  5 -node  lattice,  the  value  of  macroscopic  force  F  given  by 
Eq.  78  matches  F,  the  average  internal  force  inside  the  5-node  lattice,  as  given  by 
Eq.  28. 


To  generate  the  applied  nodal  forces,  the  relation  Eq.  58  needs  to  be  applied  to 
successive  nodes  of  the  uniformly  spaced  6-node  lattice.  A  direct  application  to 
nodes  1-3,  respectively,  gives 


Fi  = 


fi  ■  (—1)  +  /2  •  (—1)  +  /3  ■  (—1)  +  f&  ■  (—1)  +  h  •  (— 1) 


—  fi  +  H  +  h  +  fi  +  /5 


F2 


fi  '  (1)  +  fi '  (— 1)  +  h  •  (— 1)  +  fs  ■  (— 1)  +  /4  •  (— 1) 


—  /2  +  /3  +  U 


f3 


h  •  (i)  +  h  •  (i)  +  h  •  (-1)  +  h  •  (-1)  +  h  •  (-1) 


—  fs 


(79) 

(80) 

(81) 


matching  the  expressions  given  in  Eqs.  30-32  for  the  loading  of  the  uniformly 
spaced  6-node  lattice. 
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3.8  Summary  of  Uniformly  Spaced  A-Node  Lattice  Equilibrium 


The  solution  of  the  uniformly  spaced,  fully  connected  1-D-lattice  equilibrium  equa¬ 
tions  has  been  shown  using  both  force  and  energy  approaches.  Compatibility  in  the 
2  approaches  has  been  established,  whereby  it  was  shown  that  the  force  derived 
from  the  derivative  of  energy  with  respect  to  a  displacement  corresponds  to  the  av¬ 
erage  internal  force  in  the  uniformly  spaced  lattice  (and  not  to  the  net  externally 
applied  load  across  the  various  load  points). 

The  regularity  established  in  both  the  energy  and  force  equations  allows  for  the 
generalized  pattern  to  be  recognized  that  connects  the  internodal  quantities  (e^  and 
fk)  to  the  macroscopic  global  properties  (E  and  F): 

N—l 

E  =  YJ(N~^i  (82) 

i= 1 

i  n~ 1 

F  =  [  £  i(N  -  i)f,  .  (83) 

1  1=1 

The  applied  nodal  forces  and  their  net  sum,  for  a  uniformly  spaced  lattice,  can  be 
reduced  to  the  following  summations: 

N-i 

Fi  =  fm.  ,  (84) 

m—i 

N/2  N-i 

Fnet  =  ]T^/m  .  (85) 

i=  1  m=i 

A  tabulated  comparison  of  F  to  Fnet  is  given  in  the  table,  out  to  N  =  8. 

Table  A  comparison  of  global  vs.  net-applied  lattice  force  for  uniformly  spaced  lattices 


N 

F  =  F 

Fnet 

T1  -  Fnet 

2 

h 

fl 

0 

3 

fl  +  f2 

/l  +  /2 

0 

4 

+ 

COli- 

+ 

fl  +  2/2  +  h 

1/2 

5 

fl  +  |/2  +  |/3  +  fi 

fl  +  2/2  +  2/3  +  /4 

i/2  +  e>(/3) 

6 

fl  +  |/2  +  |/3  +  |/4  +  fb 

fi  +  2/2  +  3/3  +  2/4  +  F 

1/2  +  0(/3) 

7 

fi  +  I/2  +  2/3  +  2/4  +  I/5  +  fe 

fi  +  2/2  +  3/3  +  3/4  +  2/5  +  fe 

i/2 + o(/3) 

8 

fl  +  +  f/s  +  fh  +  lih  +  t/6  +  h 

fl  +  2/2  +  3/3  +  4/4  +  3/5  +  2/6  +  /7 

f/2  +  o(/3) 
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An  examination  of  the  last  column  of  the  table  shows  the  discrepancy  between  the 
global  (F)  and  net-applied  (Fnet)  lattice  force.  As  the  lattice  size  N  grows,  it  can 
be  seen  that  the  discrepancy  is  -j|/2  plus  higher  order  terms.  Considering  the  case 
where  the  lattice-spring  stiffness  is  monotonically  decreasing  beyond  a  certain  finite 
level  of  internodal  separation  (as  in  the  case  of  atomic  lattices),  one  may  conclude, 
for  large  lattices  (where  N  A>  2),  that  F  and  Fnet  converge  to  the  same  value. 

The  simply  loaded  lattices  for  IV  =  4  and  A  =  5  arc  now  examined.  However,  the 
simplifying  relations  of  Section  3.2  may  not  be  employed,  because  of  the  nonuni¬ 
form  spacing  of  the  lattice. 

3.9  The  Simply  Loaded  4-Node  Lattice 

For  this  discussion,  please  refer  to  Fig.  6.  The  simplifications  of  Section  3.2  do 
not  apply  to  this  lattice.  Nonetheless,  it  is  known  from  the  figure  that  fx  =  \x  = 
S]_a  and  lx  =  s\.  In  addition,  the  figure  informs  one  that  X/L  =  1/(3  —  2A).  A 
tabulation  of  the  nodal  energies  (Eq.  47)  gives 


2£i  =  ei_A  +  e2_A  +  £3-2A 

2F2  =  Ci  +  ei_A  +  ^2-A 

2T/3  =  €\  +  ei_A  +  C2-A 

27^4  =  ei_A  +  e2_A  +  C3_2a 
4 

E  =  'y  ^  E{  —  ei  +  2ei„A  +  2e2^A  +  C3-2A  (86) 

i=l 


A  direct  application  of  Eq.  54  reveals 


F  = 


3  -2A  L 


h  ■  1  +  2/i_a(1  -  A)  +  2/2_a(2  -  A)  +  /3_2A( 3  -  2A) 


(87) 


At  the  moment,  Eq.  87  looks  nothing  like  the  applied  lattice  force  derived  in  Sec¬ 
tion  2.7,  given  in  2  separate  forms  in  Eq.  39.  Those  2  forms  correspond  to  the 
co-equal  internal  lattice  force  at  2  separate  internal  locations,  given  by  F,  (Eq.  36) 
and  F2  (Eq.  37).  However,  consider  the  generalized  expression  for  average  internal 
lattice  force,  F,  given  in  Eq.  4.  For  the  4-node  simply  loaded  lattice  under  consid¬ 
eration,  it  becomes 


-  2Fx(l  —  A)  +  F2 
F~ - 3^2A - 


(88) 


Knowing  that  F  =  F  for  any  simply  loaded  lattice,  substituting  Eqs.  36  and  37  into 
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this  expression  produces  Eq.  87. 


This  tautological  reduction  implies  that  the  external  force  derived  from  the  energy 
approach,  Eq.  87,  corresponds  in  form  and  magnitude  to  the  average  internal  force 
derived  using  a  force  approach,  Eq.  39,  for  a  4-node,  simply  connected  1-D  lattice. 


To  generate  the  applied  boundary  force,  the  relation  Eq.  58  needs  to  be  applied  to  a 
boundary  node  of  the  4-node  simply  loaded  lattice  (see  Fig.  6).  A  direct  application 
of  Eq.  58  on  node  1  gives 


Fi  = 


fl-A  '  (  —  1)  +  /2-A  '  (  —  1)  +  /3-2A  '  (  — 1) 


(89) 


which  simplifies  to 

Fi  =  fi-A  +  /2-A  +  /3-2A  ,  (90) 

which  matches  1  of  the  expressions  given  in  Eq.  39  for  the  external  loading  of  the 
simply  loaded  4-node  lattice. 

To  generate  the  constraint  necessary  for  the  evaluation  of  A,  the  relation  Eq.  58 
needs  to  be  applied  to  an  interior  node  of  the  4-node  simply  loaded  lattice.  A  direct 
application  of  Eq.  58  on  node  2  gives 


7*2  =  0 


/l-A  •  (1)  +  A  •(-!)+ /2-A  •(-!) 


(91) 


which  simplifies  to 


fl  +  /2-A  —  /l— A  , 


which  is  the  exact  constraint  of  Eq.  38. 


(92) 


3.10  The  Simply  Loaded  5-Node  Lattice 

For  this  discussion,  please  refer  to  Fig.  7.  The  simplifications  of  Section  3.2  do  not 
apply  to  this  lattice.  Nonetheless,  it  is  known  from  the  figure  that  \x  =  \x  =  si_a 
and  \x  =  =  si.  In  addition,  the  figure  informs  one  that  X/L  =  1/(4  —  2A).  A 
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tabulation  of  the  nodal  energies  (Eq.  47)  gives 


2  El  =  ei_A+  &2-A  +  C3-A  +  ^4_2A 

2E2  =  +  Ci_A  +  62  +  C3-A 

2E/3  =  2ei  +  2e2_A 

27^4  =  Ci  +  ei_A  +  62  +  C3-A 

277.5  —  Ci_A+  e2-A  +  C3_A  +  e4_2A 

5 

E  =  'y  ^  77)  =  2ei  +  2ei_A  +  e2  +  2e2_A  +  2e3_A  +  ^4_2a 

i=  1 

A  direct  application  of  Eq.  54  reveals 


F 


1 

4  -  2A 


2fi  ■  1  +  2/1_a(1  -  A)  +  /2(2)  +  2/2_a(2  -  A) 


+  2/3_a(3  —  A)  +  /4-2a(4  —  2A) 


(93) 


(94) 


On  its  face,  Eq.  94  looks  nothing  like  the  applied  lattice  force  derived  in  Section  2.8, 
given  in  2  separate  forms  in  Eq.  45.  Those  2  forms  correspond  to  the  co-equal 
internal  lattice  force  at  2  separate  internal  locations,  given  by  F\  (Eq.  41)  and  F> 
(Eq.  42).  However,  as  before,  consider  the  average  internal  lattice  force  given  in 
Eq.  4,  for  this  5-node  lattice: 


2F2  +  2A(1-A) 

F~ - 4^2A - 


(95) 


Since  F  =  F  for  any  simply  loaded  lattice,  substituting  Eqs.  41  and  42  in  this 
expression  produces  Eq.  94. 


This  tautological  reduction  implies  that  the  external  force  derived  from  the  energy 
approach,  Eq.  94,  corresponds  in  form  and  magnitude  to  the  average  internal  force 
derived  using  a  force  approach,  Eq.  45,  for  a  5-node,  simply  connected  1-D  lattice. 


To  generate  the  applied  boundary  force,  the  relation  Eq.  58  needs  to  be  applied  to  a 
boundary  node  of  the  5-node,  simply  loaded  lattice  (see  Fig.  7).  A  direct  application 
of  Eq.  58  on  node  1  gives 


Fi  = 


fl—A  ■  (  —  1)  +  /2-A  •  (  —  1)  +  /3-A  '  (  —  1)  +  /4-2A  '  (  — 1) 


(96) 
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which  simplifies  to 


F\  —  fi-A  +  /2-A  +  /3-A  +  /4-2A  ,  (97) 


which  matches  1  of  the  expressions  given  in  Eq.  45  for  the  external  loading  of  the 
simply  loaded  4-node  lattice.  An  evaluation  of  F5  would  reveal  that  F5  =  —F\ . 


To  generate  the  constraint  necessary  for  the  evaluation  of  A,  the  relation  Eq.  58 
needs  to  be  applied  to  an  interior  node  of  the  5-node  simply  loaded  lattice.  A  direct 
application  of  Eq.  58  on  node  2  gives 


7*2  =  9 


fl—A  ‘  (1)  +  fl  ■  (“I)  +  /2  '  (“I)  +  /3-A  '  (  — 1) 


(98) 


which  simplifies  to 

fl—A  =  fl  +  f2  +  /3-A  ,  (99) 

which  is  the  exact  constraint  of  Eq.  44.  An  evaluation  of  F4  would  reveal  the  iden¬ 
tical  constraint. 


4.  Conclusion 

In  this  report,  fully  connected  1-D  lattices,  in  which  every  lattice  node  mechanically 
interacts  with  every  other  lattice  node  through  a  nonlinear  spring,  were  examined. 
No  assumptions  were  made  about  the  nature  of  the  relationship  governing  the  in- 
ternodal  spring  function  (except  by  way  of  example).  Thus,  there  is  no  presumption 
of  harmonic  lattice  behavior.  Specific  lattice  configurations  were  considered,  rang¬ 
ing  from  the  simplest  2-node  structure  to  a  6-node  lattice.  Two  constraints  were 
respectively  applied  to  each  lattice,  where  applicable:  either  the  lattice  was  con¬ 
strained  to  uniform  internodal  spacing;  or  else,  the  applied  load  was  restricted  to 
simple  loading  at  the  lattice  boundary. 

That  a  uniformly  spaced  lattice  constraint  is  incompatible  with  simple  boundary 
loading  may  not  be  surprising;  however,  nor  is  it  intuitive.  To  maintain  uniform 
spacing  in  the  lattice,  external  forces  must  be  applied  even  to  internal  nodes  of  the 
lattice.  In  contrast,  a  simple  boundary  load  upon  the  lattice  produces  a  deformation 
which  is  non-uniform.  This  occurs  because  the  aggregate  spring  constant  is  lesser 
for  nodes  near  the  lattice  boundary,  because  of  fewer  nodal  neighbors.  Reconciling 
these  discrepancies  is  essential  to  understanding  the  lattice  equilibrium  problem. 


28 


Two  methodologies  were  applied  to  each  lattice  configuration  that  was  considered: 

1)  a  force  approach,  in  which  free-body  diagrams  were  developed  that  concep¬ 
tually  cut  the  lattice  at  various  locations  in  its  internal  structure,  to  establish 
and  resolve  the  equations  of  force  equilibrium;  and 

2)  an  energy  approach,  in  which  the  potential  energy  of  every  intemodal  inter¬ 
action  was  tabulated.  With  this  approach,  both  the  global  as  well  as  the  in¬ 
dividual  nodal  force  responses  are  derived  quantities  from  the  total  potential 
energy. 

With  both  approaches,  relationships  were  established  between  the  local  internodal 
(internal)  forces  and  the  (external)  macroscopic  forces  applied  to  the  body.  The 
goal  of  the  report  was  to  confirm  the  expected,  that  these  2  approaches  yielded 
compatible  results.  That  confirmation  was  achieved  for  all  configurations  examined. 

Each  approach  has  its  strong  and  weak  aspects.  While  the  force  approach  may  be 
conceptually  easier  to  understand,  the  execution  of  it  proves  more  tedious.  In  ad¬ 
dition,  while  the  current  effort  was  limited  to  1-D  lattices,  the  force  approach  will 
need  to  further  rely  on  vector  mechanics  as  the  lattice  dimension  extends  to  2-D 
or  3-D,  thereby  adding  increasing  complexity.  The  energy  approach  does  not  suffer 
this  drawback.  The  ability  to  generalize  solutions  also  proves  much  easier  with  the 
energy  approach,  as  was  seen  in  Section  3.8. 

Key  to  achieving  the  proper  formulation  with  the  energy  approach  is  a  proper  use 
of  the  chain  rule  for  differentiation,  which  takes  the  forms  of  Eqs.  51  and  56,  as 
applied  to  the  1-D  lattice  compression  problem.  When  utilized  to  differentiate  the 
system  energy  with  respect  to  the  body  dimension,  the  resulting  global  force  equa¬ 
tion,  Eq.  54,  exhibits  a  multiplier  on  each  intemodal  force  that  is  proportional  to 
the  range  of  the  interaction.  This  particular  multiplier  is  peculiar  to  uniform  lat¬ 
tice  expansions/contractions  (as  opposed  to  distortions  from  wave  propagation),  in 
that  distal  nodes  converge  a  greater  distance  than  proximate  nodes  under  a  uniform 
decrease  in  lattice  spacing. 

An  interesting  result  of  the  study  was  the  revelation  that  the  global  force  F  that 
derives  from  the  energy  approach  was  not  necessarily  equal  to  the  net  force  Fnet 
applied  to  the  lattice.  Rather,  F  was  equivalent,  for  all  boundary  conditions  con¬ 
sidered,  to  the  average  internal  force  F  as  integrated  throughout  the  lattice.  The 
discrepancy  can  arise,  since  a  uniformly  spaced  lattice  (for  N  >  3)  receives  applied 
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loads  to  internal  nodes,  which  means  that  such  loads  do  not  deform  the  complete 
lattice  structure,  but  only  parts  thereof.  However,  as  shown  in  section  3.8,  the  dis¬ 
crepancy  will  tend  to  0  as  the  lattice  size  N  grows,  assuming  the  intemodal  lattice 
stiffness  decreases  with  increasing  nodal  separation,  as  in  the  case  of  atomic  lat¬ 
tices. 


30 


5.  References 


1.  Harmon  PM.  Energy,  force,  and  matter:  The  conceptual  development  of 
nineteenth-century  physics.  Cambridge:  Cambridge  University  Press;  2005. 

2.  Zienkiewicz  OC,  Taylor  RL.  The  finite  element  method.  5th  ed.  Oxford: 
Butterworth-Heinemann;  2000.  First  published  by  McGraw-Hill  in  1967. 

3.  Bathe  KJ.  Finite  element  procedures  in  engineering  analysis.  Englewood  Cliffs: 
Prentice  Hall;  1982. 

4.  Sears  FW,  Zemansky  MW,  Young  HD.  University  physics.  Fifth  ed.  Reading: 
Addison- Wesley;  1977. 

5.  Beer  FP,  Johnston  Jr.  ER.  Vector  mechanics  for  engineers:  statics.  Third 
ed.  New  York:  McGraw  Hill;  1977. 

6.  Segletes  SB.  The  vibrational  stiffness  of  an  atomic  lattice.  Aberdeen  Proving 
Ground  (MD):  Army  Research  Faboratory  (US);  1998  Sep.  Report  No.:  ARF- 
TR-1757. 

7.  Segletes  SB.  Vibrational  stiffness  of  an  atomic  lattice.  In:  Furnish  MD,  Chhabil- 
das  FC,  Hixson  RS,  editors.  Shock  Compression  of  Condensed  Matter — 1999; 
Vol.  1;  1999  Jun27  -  Jul  2;  Snowbird,  UT.  p.  161-164.  AIP  Conference  Proceed¬ 
ings  505. 

8.  Segletes  SB.  Classical  methods  for  frequency-based  equations  of  state.  Ab¬ 
erdeen  Proving  Ground  (MD):  Army  Research  Faboratory  (US);  2007  Mar.  Re¬ 
port  No.:  ARF-RP-166.  (a  reprint  from  Recent  Res.  Devel.  Phys.  Chem.  Solids, 
vol  1.  pp.  1-26,  2002). 


31 


Intentionally  left  blank. 


32 


List  of  Symbols,  Abbreviations,  and  Acronyms 


A  -  the  nominal  lattice  spacing  under  the  given  load  ( i.e .,  the  actual  lattice  spac¬ 
ing  “far”  from  the  lattice  boundary). 

A  -  the  relative  amount  by  which  a  given  lattice  cell  length  may  deviate  from  the 
nominal  value  of  A,  owing  to  a  change  in  lattice  stiffness  near  the  boundary 
of  the  lattice.  Its  value  is  limited  to  real  numbers  in  the  range  A  <  1. 

€i  -  the  infinitesimal  fraction  of  A  by  which  a  lattice  node  is  moved  in  order  to 
assess  a  change  in  the  net  potential  energy  of  the  lattice. 

N  -  the  total  number  of  nodes  in  the  1-D  lattice. 

L  -  the  total  length  of  the  macroscopic  lattice,  L  ps  (N  —  1)A.  However,  when 
the  lattice  is  constrained  to  uniform  spacing,  the  approximation  is  exact. 

Sk  -  shorthand  for  internodal  distance  in  wavelengths,  such  that  Sk  =  kX.  The 
index  k  may  take  on  non-integer  values  if  the  lattice  spacing  varies  from  the 
uniformly  spaced  configuration. 

™x-  the  (positive)  physical  separation  distance  between  node  pair  i^rn  of  the 
lattice.  For  the  special  case  of  a  uniformly  spaced  lattice,  ™x  =  S|m_j|  = 
\m  —  i\X.  Whereas  the  notation  "'x  provides  the  specificity  of  both  nodes  in 
the  node  pair,  the  notation  Sk  allows  the  convenience  of  tracking  only  the 
separation  distance  in  node  pair  iv-rn,  without  the  requirement  to  know  the 
values  of  i  and  m,  specifically. 

fk  -  the  internal  spring  force,  using  the  shorthand  that  fk  =  f(sk )  =  f(kX)  (note 
that  /o  =  0).  The  index  k  may  take  on  non-integer  values. 

Fj  -  the  external  Cartesian  load  applied  to  node  i.  Positive  denotes  a  force  in 
the  rightward  direction,  while  negative  denotes  a  force  to  the  leftward  di¬ 
rection.  Note  that,  because  of  the  symmetric  equilibrium  being  imposed, 

Fj  =  —FN+1_j. 

Fnet  -  the  net  externally  applied  load,  Fnet  =  Fj- 

F  -  when  subject  to  simple  boundary-loading,  the  internal  force  at  any  point  in 
the  lattice. 
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Fj  -  the  internal  force  in  a  lattice  (positive  implies  compression)  that  exists  be¬ 
tween  nodes  j  and  (j  +  1)  (e.g.,  F\  is  the  internal  force  between  the  node 
pair  1tA2).  Note  that,  because  of  the  symmetric  equilibrium  being  imposed, 
Fj  =  FN_j.  When  a  lattice  is  subject  to  simple  boundary  loading,  F}  =  F. 

F  -  the  average  internal  load  in  the  lattice,  given  generally  by 

F  —  y  f  F(x)  dx 
L  Jo 

TV- 1  IN- 1 

=  E  b  S+F  /  E  i+F  ■ 

j= 1  /  3=1 


In  the  case  of  a  uniformly  spaced  lattice,  where  j+jx  =  A,  the  above  reduces 
to 


F 


1 


N-  1 


TV— 1 


F  -  when  using  the  energy  method,  F  is  the  derivative  of  the  system  mechanical 
energy  with  respect  to  a  change  in  the  body  dimension.  It  was  determined 
that,  in  the  case  of  both  uniformly  spaced  as  well  as  simply  loaded  lattices, 
F  is  equivalent  to  the  average  internal  force,  F. 


ek  -  an  atom’s  potential  energy  arising  from  an  interaction  with  an  atom  at  a 
distance  of  Sk ■  Thus,  =  e(k\)  Note  that  e0  =  0.  The  index  k  may  take  on 
non-integer  values. 


Et  -  Net  energy  potential  associated  with  atom  i  in  the  lattice. 

E  -  the  total  energy  potential  of  the  lattice,  E  =  J2iLi  - 
Ti  aAz2  -  the  node  pair  i ] ,  i2  through  which  an  internal  nonlinear  spring  force  acts. 
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